Zienolddiny et al. Breast Cancer Research 2013, 15:R53 
http://breast-cancer-research.eom/content/15/4/R53 



Breast Cancer 

I ?1 RESEARCH 



RESEARCH ARTICLE 



Open Access 



Analysis of polymorphisms in the circadian- 
related genes and breast cancer risk in 
Norwegian nurses working night shifts 

Shanbeh Zienolddiny 1 *, Aage Haugen 1 , Jenny-Anne Sigstad Lie 2 , Helge Kjuus 2 , Kristine Haugen Anmarkrud 1 and 
Kristina Kjaerheim 3 



Abstract 

Introduction: Some studies have suggested that night work may be associated with an increased risk of breast 
cancer in nurses. We aimed to explore the role of circadian gene polymorphisms in the susceptibility to night 
work-related breast cancer risk. 

Methods: We conducted a nested case-control study of Norwegian nurses comprising 563 breast cancer cases and 
619 controls within a cohort of 49,402 Norwegian nurses ages 35 to 74 years. We studied 60 single-nucleotide 
polymorphisms (SNPs) in 17 genes involved in the regulation of the circadian rhythm in cases and controls. The 
data were analyzed in relation to the two exposure variables "maximum number of consecutive night shifts ever 
worked" and "maximum number of consecutive night shifts worked for at least 5 years." The odds of breast cancer 
associated with each SNP was calculated in the main effects analysis and in relation to night shift work. The 
statistically significant odds ratios were tested for noteworthiness using two Bayesian tests: false positive report 
probability (FPRP) and Bayesian false discovery probability (BFDP). 

Results: In the main effects analysis, CC carriers of rs4238989 and GG carriers of rs3760138 in the AANAT gene had 
increased risk of breast cancer, whereas TT carriers of BMAL1 rs2278749 and TT carriers of CLOCK rs3749474 had 
reduced risk. The associations were found to be noteworthy using both the FPRP and BFDP tests. With regard to 
the effect of polymorphisms and night work, several significant associations were observed. After applying FPRP 
and BFDP in women with at least four night shifts, an increased risk of breast cancer was associated with variant 
alleles of SNPs in the genes AANAT (rs3760138, rs4238989), BMAL1 (rs2290035, rs2278749, rs969485) and ROR-b 
(rs3750420). In women with three consecutive night shifts, a reduced risk of breast cancer was associated with 
carriage of variant alleles of SNPs in CLOCK (rs3749474), 8MAL1 (rs2278749), BMAL2 (rs2306074), CSNK1E (rs5757037), 
NPAS2 (rsl 7024926), ROR-b (rs3903529, rs3750420), MTNR1A (rs1 31 1 1 3549) and PER3 (rs1012477). 

Conclusions: Significant and noteworthy associations between several polymorphisms in circadian genes, night 
work and breast cancer risk were found among nurses who had worked at least three consecutive night shifts. 
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Introduction 

Breast cancer in women is the second leading cause of 
mortality [1]. The etiology of breast cancer involves several 
risk factors, such as early age at menarche, late age at 
menopause, late age at first birth, nulliparity, hormonal 
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factors, ionizing radiation exposure and alcohol consump- 
tion [2,3]. However, these known risk factors do not fully 
explain the high incidence of breast cancer. Genetic fac- 
tors are also involved in the etiology of breast cancer, 
including rare high-risk mutations (BRCA1 and BRCA2), 
more moderate susceptibility variants (CHEK2, ATM) and 
several as yet unidentified common susceptibility variants 
associated with low to moderate increased risk [4,5] . 
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At present, there is not sufficient evidence for any 
occupational exposure to be classified as carcinogenic in 
relation to breast cancer. However, night work is a 
potential risk factor that has been found in many studies 
to be associated with an elevated risk of breast cancer 
[6-12]. A recent study reported that 4.6% of all female 
breast cancers could be attributed to occupational expo- 
sure to shift work [13]. This is important, given the 
large proportion of women engaged in shift work and 
the very high incidence of female breast cancer in gen- 
eral. Although the majority of epidemiological studies 
are suggestive of a positive association, not all studies 
show a consistent association between night work and 
increased risk of breast cancer [14-16]. 

Night work was classified as a Group 2A carcinogen by 
the International Agency for Research on Cancer (IARC) 
in 2007 [17,18]. The hypothesis of circadian disruption 
by light at night (LAN) has received particular attention 
in this evaluation [19]. Disruption of the circadian 
rhythm has been hypothesized to increase cancer risk 
through its impact on the production of the circadian 
hormone melatonin [20,21]. A potential mechanism 
could be through suppression of the nocturnal produc- 
tion of melatonin, caused by exposure to LAN, which 
may also influence the patterns of sex hormone produc- 
tion and in turn influence breast cancer risk [22,23]. The 
levels of melatonin are regulated through its biosynthesis 
from the amino acid tryptophan, which is mediated by 
the arylalkylamine 7V-acetyltransferase (AANAT) protein 
or melatonin receptors MTNR1A (MT1) and MTNR1B 
(MT2) and possibly the retinoic acid-related orphan 
receptors ROR-a and ROR-b [24]. An antiproliferative 
effect of melatonin has also been shown in human breast 
cancer cell lines. In estrogen receptor-positive (ER + ) 
human breast cancer cells, melatonin suppresses both ER 
mRNA expression and estrogen-induced transcriptional 
activity of the ER [25,26]. Melatonin also regulates the 
transcriptional activity of other enzymes involved in 
estrogen metabolism and the expression of circadian and 
circadian-targeted genes. Recently, melatonin was found 
to induce apoptosis [27] and to reduce the levels of oxi- 
dative stress [28,29]. Other mechanisms, such as impact 
of sleep disruption on human physiological processes, 
may also be involved [30], although not all the results of 
studies on the impact of sleep disruption on breast can- 
cer risk have been consistent [31,32]. 

The circadian clock is regulated by several proteins, such 
as aryl hydrocarbon receptor nuclear translocator-likes 
(ARNTL or BMAL1, ARNTL2 or BMAL2) also known as 
Mop3, circadian locomotor output cycle kaput (CLOCK), 
neuronal PAS domain protein 2 (NPAS2), periods (PERI, 
PER2 and PER3), cryptochromes (CRY1 and CRY2), 
TIMELESS, ROR-a and ROR-b, that interact with each 
other in various transcriptional and translational feedback 



loops [33]. The main signal is the daily light-dark rotation, 
which is mediated by light-sensing opsin proteins such as 
opsin 4 (OPN4) [34]. 

The circadian genes may affect cancer susceptibility 
through effects on biological pathways that regulate 
DNA damage and repair, carcinogen metabolism and/or 
detoxification, cell growth and cell death [35]. It has 
been shown that expression of the circadian gene PERI 
may be downregulated in tumors and tumor cell lines 
compared to normal cells [36]. The PER2 gene has been 
shown to induce apoptosis through upregulation of the 
tumor suppressor gene TP53 and the proapoptotic gene 
Bax [37]. A PER2 gene mutation (S662G) also leads to 
enhanced resistance to X-ray-induced apoptosis and 
increased RAS-mediated oncogenic transformation, and 
it affects tumorigenesis in mice [38]. 

To date, there are no genetic susceptibility biomarkers 
identifying women who are particularly susceptible to 
breast cancer related to exposure to night work. In parti- 
cular, knowledge of how night work, together with com- 
mon genetic variants in genes regulating the circadian 
rhythm, may affect breast cancer risk is limited [39] . The 
aim of this study was to analyze a number of polymorph- 
isms in the core circadian genes and genes involved in 
melatonin biosynthesis and melatonin bioavailability as 
potential biomarkers of breast cancer risk in night 
workers. 

Methods 

Case patients and controls 

The details of study design, data collection, recruitment 
of cases and controls have been described previously by 
Lie et al. [8,9]. Briefly (see Additional File 1, Figure SI), 
a nested case-control study of 1,182 women ages 35 to 
74 years was performed within an established national 
cohort of 49,402 Norwegian nurses who graduated 
between 1914 and 1985. All of the women had been 
diagnosed with histologically confirmed invasive breast 
cancer between 1990 and 2007 (n = 1,132), were alive as 
of February 2009 (n = 943), consented to be interviewed 
(« = 706) and had an occupational history as a nurse of 
at least 1 year (« = 699). The controls were frequency- 
matched to cases by year of diagnosis and in 5-year age 
groups. They were all cancer-free at and prior to the year 
of diagnosis of the case, were alive as of February 2009 
(« = 1,384), had consented to an interview [n = 900) and 
worked as nurses for minimum of 1 year (n = 895). 
Breast cancer cases were identified by linkage to the 
Cancer Registry of Norway by personal identification 
number given to all Norwegians alive since 1960. The 
Cancer Registry of Norway has virtually complete records 
of all cancer cases in the Norwegian population since 
January 1953 [40]. Both cases and controls gave their full 
informed consent that their information could be used 
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and published for research purposes, given that their per- 
sonal details would remain anonymous. The study was 
approved by the Regional Committee for Medical and 
Health Research Ethics, South-East region (S-08430a, 
2008/10453). 

Assessment of night work 

Full details on the exposure assessment were described 
in a recent study [9]. Briefly, different exposure metrics 
of night work were computed on the basis of work his- 
tory. The term night work in this study includes working 
periods from rotating as well as permanent night sche- 
dules. A "night shift" was defined as a shift including 
work between 12 a.m. and 6 a.m. For each job, we 
assessed information gathered for the year of start and 
end, workplace (hospital, other institution and/or 
others), proportion of full job and work schedule (only 
day shifts, only night shifts or both day and night shifts). 
If the job included night shifts, we also asked about the 
number of night shifts per month, and number of con- 
secutive night shifts. Two surrogate measures for night 
work were applied: maximum number of consecutive 
night shifts in all jobs (<3, 3 and >4) and maximum 
number of consecutive night shifts the nurse had 
worked during a minimum of 5 years (<3, 3 or >4). 

We focused our analysis on the variable consecutive 
night shift because, in our previous analysis [9], the con- 
secutive night shift variable that takes duration into 
account was the only variable found to be significantly 
associated with an increased risk of breast cancer. In 
that previous analysis, a positive trend was seen with an 
increasing number of consecutive nights, albeit that a 
significant association was observed only when compar- 
ing those with at least 5 years of work in a job including 
six or more consecutive night shifts with those who 
never worked nights. 

Collection of saliva and DNA extraction 

DNA from 563 saliva samples taken from cases and 619 
samples obtained from controls were isolated using the 
Oragene-DNA isolation kit (DNA Genotek Inc, Kanata, 
ON, Canada) as instructed by the manufacturer, with 
some modifications. Briefly, 1 ml of saliva was trans- 
ferred to an Eppendorf tube, and reagents provided by 
the supplier were added. DNA was successfully 
extracted from all saliva samples that were received. 

Selection of polymorphisms and genotyping 

The single-nucleotide polymorphism (SNP) selection strat- 
egy was based primarily on a candidate gene approach 
using information from a search of relevant studies pub- 
lished by 2009. We searched the literature for polymorph- 
isms in the core circadian genes and melatonin biosynthesis 
and/or bioavailability that had been investigated for an 



association with cancer risk. Since there were few poly- 
morphisms that had been investigated by 2009, we also 
included polymorphisms using one or more of the follow- 
ing criteria: assumed functionality (located in the regulatory 
regions, for example, 3'-UTR, 5'-UTR or amino acid 
change), cancer genetic marker of susceptibility (CGEMS) 
for breast cancer in epidemiological studies (candidate gene 
or genome-wide association studies (GWASs)), candidate 
tag SNP with r 2 > 0.80 and minor allele frequency (MAF) > 
5%, and location in the relevant circadian, melatonin and 
light signaling pathways. Because of the constraint of the 
project budget, and hence limitations on the number of 
polymorphisms that could be placed on the chip, 60 SNPs 
were successfully designed and genotyped (see Additional 
File 2, Tables SI and S2). The genotyping success rate 
was similar in cases and controls: at least 0.80 in cases and 
0.82 in controls. All SNPs except for one SNP in the CRY2 
gene (rsl2364060) and one SNP in the MTNR1A gene 
(rs2 119882) were in Hardy- Weinberg equilibrium. Two 
previously reported SNPs, rs35670208 in the ARNTL2/ 
BMAL2 gene and rs28936679 in the AANAT gene, respec- 
tively, were found to be nonpolymorphic in this population. 
Genotyping was performed using the iPLEX genotyping 
assay (Sequenom Inc, San Diego, CA, USA) at the Centre 
for Integrative Genetics (CIGENE) genotyping core facility 
at the University for Environment and Biosciences, As, 
Norway. 

Statistical methods 

Pearson's j 2 test was used to assess case-control differ- 
ences in the frequency distribution of categorical variables, 
and a nonparametric Mann-Whitney test was used to 
assess differences in median age between cases and con- 
trols. On the basis of existing knowledge of factors asso- 
ciated with breast cancer, the following factors were 
regarded as potential confounders: age at time of diagnosis 
(as a continuous variable in 5 -years age groups), period of 
diagnosis (1996 to 1999 and 2000 to 2007), parity (0, 1 or 
2 or >3 children), family history of breast cancer in mother 
and/or sister (no or yes), use of hormone treatment during 
the 2 years before the reference year (no or premenopau- 
sal, or yes) and frequency of alcohol consumption at time 
of diagnosis (less than twice weekly or at least twice per 
week). These variables were included in the final logistic 
regression model to obtain adjusted odds ratios (ORs) for 
developing breast cancer associated with the frequency of 
polymorphisms alone (main effects analysis), the poly- 
morphisms and night work (stratified by number of conse- 
cutive night shifts). On the basis of the results of our 
previous study showing associations between intensity of 
night work and breast cancer risk, we examined the simul- 
taneous effects of polymorphisms and night work by com- 
paring the frequency of genotypes of each polymorphism 
in cases and controls who had worked either at least four 
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and three consecutive night shifts with a reference group 
who had worked only day shifts or a maximum of two 
consecutive night shifts. We compared the frequency of 
SNPs in subgroups defined by both intensity (number 
of consecutive nights) and duration (number of years) of 
night work. In our previous study [9], the consecutive 
night shift variable that takes duration into account was 
the only variable found to be significantly associated with 
an increased risk of breast cancer, that is, a significant 
association was observed only when individuals with at 
least 5 years of work experience in a job including six or 
more consecutive night shifts were compared with those 
who had never worked nights. A positive trend was seen 
with increasing number of consecutive nights, however, 
thus, to avoid small numbers of subjects in the stratified 
analyses in the present study, women who had worked at 
least four consecutive night shifts were classified as the 
most exposed group. The relative risk of breast cancer 
associated with each genotype was estimated as OR with 
95% confidence interval (CI) using unconditional logistic 
regression models. For determination of interactions 
between the shift work and SNP genotypes, the product of 
shift work x genotype variables was added to the model. 
Interaction P values {P interaction) are reported for the 
joint effects of shift work and genotype. All reported 
P values are two-sided. P < 0.05 was considered statisti- 
cally significant. All statistical analyses were carried out 
using SPSS release 21 software (SPSSInc, Chicago, IL, 
USA). 

Multiple testing and false discovery rates 

To evaluate the possibility of false-positives due to mul- 
tiple testing, we applied two Bayesian statistical tests, 
false positive report probability (FPRP) and the Bayesian 
false discovery probability (BFDP), that have been used 
frequently to test noteworthiness of findings in genetic 
association studies [41-43]. The FPRP statistical tool was 
introduced by Wacholder et al. [42] as a means to 
assess whether the strength of an association is note- 
worthy in the case of multiple hypothesis testing. In 
addition, we used the BFDP, which is a recent further 
development of FPRP [41]. We performed both FPRP 
and BFDP tests for each SNP. For both tests, we chose 
the noteworthiness thresholds (FPRP 0.2 and BFDP 0.8) 
following the original papers [41,42]. We considered 
noteworthy only those associations that passed both 
tests. We have set low prior probability values compared 
to what commonly has been used [44]. 

Results 

The study design is a nested case-control study (see 
Additional File 1, Figure SI). The original cohort study 
of Norwegian nurses was described in detail by Lie et al. 
[8,9]. The demographic characteristics of cases and 



controls (Table 1) are not significantly different, except 
for the presence of breast cancer in mother or sister. 
The exposure parameters and the number of cases and 
controls in each exposure group are shown in Table 2. 
The selected circadian genes, chromosome locations and 
the protein function of each respective gene are shown 
in Additional File 2, Table SI. The SNPs with corre- 
sponding reference SNP (rs) numbers, SNP types, base 
changes, MAF, Hardy- Weinberg-derived P values and 
genotyping rates are given in Additional File 2, Table S2. 

Sixty SNPs in 17 genes (see Additional File 2, Tables 
SI and Table S2) from the circadian signaling pathway 
were analyzed for association with breast cancer risk. As 
shown in Additional File 2, Table S2, the MAFs of the 
SNPs in controls were similar to those reported pre- 
viously in Caucasian populations, indicating the validity 
of the genotyping methodology [45-47]. 

The analysis of main effects, without regard to night 
work (significant associations are shown in Table 3), 
showed that four SNPs were associated with breast can- 
cer risk: CC carriers of the AANAT rs4238989 with an 
OR of 1.43 (95% CI, 1.01- 2.01) and GG carriers of 
AANAT rs3760138 with an OR of 1.42 (95% CI, 1.01- 
2.01) had increased risk of breast cancer. The carriers of 
variant TT genotype of SNP rs2278749 in the BMAL1 
gene with an OR of 0.45 (95% CI, 0.23- 0.90) and the 
homozygote genotype TT of rs3749474 in the CLOCK 
gene with an OR of 0.64 (95% CI, 0.45- 0.92) had 
reduced risk of breast cancer. All of these associations 
were found to be noteworthy for both FPRP (<0.2) and 
BFDP (<0.8). 

Compared to the reference group (Table 4), in the high- 
est exposed group with at least four consecutive night 
shifts, CLOCK rsl 1133373 (OR 2.75, 95% CI 1.07- 7.09), 
ROR-b rs3903529 (OR 1.87, 95% CI 1.0- 3.61), two SNPs 
in the AANAT gene (rs3760138 (OR 1.82, 95% CI 1.13- 
2.93) and rs4238989 (OR 1.62, 95% CI 1.0- 2.68)) and 
Af77v7?i£-rsl0830963 (OR 2.48, 95% CI 1.06- 5.8) were 
associated with increased risk of breast cancer. In the 
same group, two SNPs in the CRY2 gene, rsll038689 (OR 
0.42, 95% CI 0.17- 1.01) and rsl401417 (OR 0.31, 95% CI 
0.10- 0.94), were associated with reduced risk. In the 
group with three night shifts, a reduced risk of breast can- 
cer was associated with the presence of one variant allele 
(heterozygote) or two variant alleles (homozygote) for 
SNPs in CLOCK rslll33373 (OR 0.50, 95% CI 0.25- 0.99), 
CLOCK rs3749474 (OR 0.56, 95% CI 0.31- 0.99); ROR-b 
rs3903529 (OR 0.57, 95% CI 0.34- 0.95), ROR-b rs7022435 
(OR 0.16, 95% CI 0.05- 0.55); ROR-b rs7022435 (OR 0.17, 
95% CI 0.03- 0.84); ROR-b rs3750420 (OR 0.40, 95% CI 
0.18- 0.88) and MTNR1A rs 13113549 (OR 0.50, 95% CI 
0.28- 0.90). 

The cases and controls were then grouped according 
to both intensity and duration of night work, that is, 
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Table 1 Characteristics of controls and breast cancer cases. 


Parameter 


Cases (n = 563) 


Controls (n = 619) 


P value 


Aae (vrf 






0.5 1 b 


Mean (SD) 


52.40 (7.88) 


52.37 (8.18) 




f v itrw i q i i 


50 


55 




Min/max 


34/70 


34/70 




Selection period 








1 996 to 1 999/2000 to 2007 (n) 


1 95/368 


221/398 


0.72 c 


Breast cancer in first-degree family d (n) 


557 


615 




Yes/no 


1 04/453 


54/561 


<0.001 


Hormone therapy in past 2 yr e (n) 


552 


605 




Yes/no 


1 27/425 


121/484 


0.22 c 


Alcohol use f 


n = 563 


n = 619 




Yes/no 


43/520 


37/582 


0.26 c 


Children (n) 


563 


619 




0/>1 


77/486 


81/538 


0.80 c 


a Age at selection for controls and age at diagnosis for case? 


. b Derived from Mann-Whitney tMest for two independent samples 


(two-sided) assuming the 


medians are the same across cases and controls. c Pearson's 


X test (two-sided). 


Cancer in sister or mother. e Hormone therapy 


among postmenopausal women 


only. f Alcohol consumption at least twice weekly. 









maximum number of consecutive nights they had ever 
worked for at least 5 yr. In subjects with at least four 
consecutive night shifts, increased risk of breast cancer 
was associated with variant alleles of BMAL1 (rs2290035 
(OR 1.91, 95% CI 1.08- 3.37; and OR 1.95, 95% CI 0.97- 
3.90) and rs969485 (OR 1.64, 95% CI 1.03- 2.61; and 
OR 1.63, 95% CI 0.96- 2.76)) and ROR-b (rs3750420 
(OR 1.60, 95% CI 1.0- 2.56), rs3903529 (OR 2.77, 95% 
CI 0.98- 7.86)) and MTNR1B rsl0830963 (OR 9.12, 95% 
CI 1.14- 73.09). In the group with three consecutive 
night shifts, the SNP rsll695472 in the PER2 gene (OR 
2.69, 95% CI 1.08- 6.73) was associated with increased 
risk of breast cancer. The SNPs BMAL1 rs2278749 (OR 
0.51, 95% CI 0.29- 0.87), BMAL2 rs2306074 (OR 0.30, 
95% CI 0.17- 0.75), CLOCK rs3749474 (OR 0.24, 95% CI 
0.05 to 1.18), CSNK1E (rsl534891 (OR 0.55, 95% CI 
0.30- 1.01), rs5757037 (OR 0.67, 95% CI 0.45- 0.99)), 



NPAS2 rsl7024926 (OR 0.33, 95% CI 0.13- 0.84; and 
OR 0.24, 95% CI 0.09- 0.67); PER3 rsl012477 (OR 0.42, 
0.24- 0.72; and OR 0.38, 95% CI 0.18- 0.66), ROR-b 
(rs7022435 (OR 0.25, 95% CI 0.07- 0.89), rs3750420 (OR 
0.37, 95% CI 0.14- 0.96; and OR 0.33, 95% CI 0.11- 1.0)) 
and MTNR1A rsl31113549 (OR 0.47, 95% CI 0.25- 
0.91) were associated with decreased breast cancer risk 
(Table 5). 

After applying FPRP and BFDP tests, in subjects with 
four or more night shifts (Table 4), increased risk of breast 
cancer associated with homozygote carriers of the variant 
alleles of SNPs rs3760138 and rs4238989 in the AANAT 
gene was found to be noteworthy. In subjects with three 
consecutive night shifts, the odds of reduced risk of breast 
cancer associated with carriage of variant alleles of SNPs 
in CLOCK (rs3749474), ROR-b (rs3903529, rs3750420) 
and MTNR1A (rsl31113549) were found noteworthy. 



Table 2 Occupational exposure data and number of cases and controls in each exposure group. 



Maximum number of consecutive night shifts in all jobs that nurse has worked 3 No. of No. of 

cases b controls b 



Only day work or one or two consecutive night shifts (reference) 


1 70 (227) 


186 (275) 


Three consecutive night shifts 


156 (181) 


189 (272) 


Four or more consecutive night shifts 


239 (291) 


244 (348) 


Total 


563 (699) 


619 (895) 


Intensity of night work: Maximum number of consecutive night shifts the nurse has ever worked for 5 or 






more years' 






Only day work or one or two consecutive nights for 5 or more years (reference) 


326 (421) 


362 (524) 


Three consecutive night shifts for 5 or more years 


125 (147) 


154 (227) 


Four or more consecutive night shifts for 5 or more years 


112 (131) 


103 (144) 


Total 


563 (699) 


619 (895) 



a Maximum number of consecutive night shifts in all jobs: 0 to 2, 3 or >4. b Number of nurses who had consented to send saliva for genetic analysis (numbers in 
parentheses show the total number of cases and controls interviewed in each exposure group). c Maximum number of consecutive night shifts the nurse has ever 
worked for 5 or more years. 
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Table 3 Odds ratios and 95% confidence intervals for polymorphisms with statistically significant associations with 
risk of breast cancer. 



Genes and SNPs 



Frequency of genotypes 3 



OR (95% Cl) b F 



FPRP" BFDP 



BMAL1 

rs2278749 
CC 
CT 
TT 

CLOCK 

rs3749474 

CC 

CT 

TT 

AANAT 

rs3760138 

TT 

TG 

GG 

rs4238989 
GG 
GC 
CC 



No. of cases 

375 
141 
13 



223 
251 
68 



149 

269 
128 



161 
263 
120 



No. of controls 

381 
164 
30 



228 
260 



315 
106 



20:! 



10--! 



1.0 

0.93 (0.70- 1 .22) 

0.45 (0.23- 0.90) 0.02 

Ptrend = 0.07 



1.0 

1.01 (0.78- 1.31) 

0.64 (0.45- 0.92) 0.02 

P trend = 0.03 



1.0 

1 .00 (0.75- 1 .32) 

1.42 (1.01- 2.01) 0.05 

Ptrend = 0.07 

1.0 

1.13 (0.86- 1.49) 

1.43 (1.01- 2.01) 0.04 

Ptrend = 0.13 



0.16 0.21 



0.14 0.11 



0.13 0.18 



0.11 0.16 



a Number of cases and controls with successful genotyping data. b Odds ratios (ORs) derived from unconditional logistic regression and adjusted for age at 
diagnosis (cases) or day of reference (controls), selection period, alcohol use, hormone replacement therapy for the past 2 years, number of children and history 
of breast cancer in first-degree family (mother or sister). c Pearson's i 2 test (two-sided) comparing frequency of the homozygote variant genotypes with 
homozygote common genotypes (reference). d The false-positive report probability (FPRP) was calculated as described by Wacholder et al. [42]. FPRP < 0.2 was 
considered noteworthy. e The Bayesian false-discovery probability (BFDP) was calculated as described by Wakefield ef al. [41]. BFDP < 0.8 was considered 
noteworthy. For more details on FPRP and BFDP, see Methods section. P t rend values for the Cochran-Armitage trend test were calculated from 2xC tables using 
frequency of cases and controls. 



With respect to both intensity and duration of night shift 
work for at least 5 yr (Table 5), in subjects with at least 
four consecutive night shifts, increased risk of breast cancer 
was associated with variant alleles of BMAL1 (rs2290035, 
rs969485) and ROR-b (rs3750420). In the subjects with 
three night shifts, SNPs in BMAL1 (rs2278749), BMAL2 
(rs2306074), CSNK1E (rs5757037), NPAS2 (rsl7024926), 
PER3 (rsl012477) and MTNR1A (rsl31113549) were asso- 
ciated with decreased breast cancer risk. Haplotype analysis 
of the polymorphisms localized on the same chromosome 
did not result in any significant associations (data not 
shown). 

Discussion 

In the main effects analysis, two SNPs in BMAL1 and 
CLOCK were associated with risk of breast cancer. There 
are no functional data to explain how each SNP in each of 
the genes may affect the risk for breast cancer. SNPs may 
change transcription binding affinity of the BMAL1/ 
CLOCK complex that is necessary for expression of the 
target genes. For example, a reduced affinity of BMAL1/ 



CLOCK to CLOCK promoter due to polymorphisms in 
these genes may lead to lower expression of CLOCK pro- 
teins. A reduction in the levels of another core circadian 
protein (CRY2) in breast epithelial cells in culture has 
been shown to lead to higher accumulation of DNA 
damage after exposure to chemical mutagens [48,49]. 
Downregulation of CLOCK has also been reported to 
cause alterations in expression of genes involved in breast 
carcinogenesis [50]. 

Several SNPs in the core circadian genes had noteworthy 
associations with breast cancer risk. Women with the 
highest number of successive night shifts and carrying at 
least one variant allele of SNPs in the two core circadian 
genes BMAL1 (rs2290035, rs969485) and ROR-b 
(rs3903529, rs3750420) had a noteworthy increased risk of 
breast cancer. Although not noteworthy, there was also a 
2.75 increased risk of breast cancer in women who had 
worked four or more night shifts and carrying two variant 
G alleles of CLOCK rsllll33373 compared to women 
with the same genotype in the reference group. Interest- 
ingly, women with three successive night shifts and 
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Table 4 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each 
polymorphism among women with less than two consecutive night shifts compared to women with three consecutive 
night shifts or women with four or more consecutive night shifts. 

Genes and SNPs Number of consecutive night shifts that nurse has worked in all jobs 

0 to 2 consecutive night Three consecutive night shifts >4 consecutive night shifts 

shifts 

FPRP d BFDP e 





Ca a /Co a 


OR (95% Cl) b 


Ca a /Co a 


OR (95% Cl) b P inter c 


FPRP d BFDP e 


Ca a /Co a 


OR (95% Cl) b P lnu 


CLOCK 
















rsl 1133373 
















CC 


72/77 


1.0 (reference) 


71/68 


0.82 (0.52 - 1 .30) 




1 03/99 


1.29 (0.79 - 2.11) 


CG 


76/72 


0.98 (0.61 -1.57) 


58/82 


0.77 (0.53 - 1.13) 




93/114 


0.87 (0.63 - 1 .20) 


GG 


14/25 


1 .05 (0.49 - 2.29) 


13/29 


0.50 (0.25 - 0.99) 0.05 
P trend f = 0.04 


0.22 0.23 


28/23 


1.17 (0.65 - 2.10) 
P trend f = 0.95 


CC 


72/77 


1.0 (reference) 


71/68 


0.82 (0.52 - 1 .30) 




1 03/99 


1.29 (0.79 - 2.11) 


CG 


76/72 


1.0 (reference) 


58/82 


0.88 (0.56 - 1 .39) 




93/114 


0.79 (0.48 - 1 .30) 


GG 


14/25 


1.0 (reference) 


13/29 


0.83 (0.34 - 2.04) 




28/23 


2.75 (1.07 - 7.09) 


rs3749474 
















CC 


67/71 


1.0 (reference) 


60/66 


0.95 (0.57 - 1.57) 




96/91 


1.14 (0.72 - 1.80) 


CT 


77/76 


1.12 (0.69-1.81) 


73/78 


1.01 (0.70 - 1.46) 




101/106 


1 .03 (0.74 - 1 .42) 


tt 


16/30 


0.55 (0.27-1.14) 


20/37 


0.56 (0.31 - 0.99) 0.05 

Ptrend' = 0.30 


0.18 0.21 


32/41 


0.85 (0.52 - 1 .39) 
P trend f = 0.54 


CC 


67/71 


1.0 (reference) 


60/66 


0.95 (0.57 - 1.57) 




96/91 


1.14 (0.72 - 1.80) 


CT 


77/76 


1.0 (reference) 


73/78 


0.93 (0.58 - 1 .48) 




101/106 


0.98 (0.63 - 1.51) 


TT 


16/30 


1.0 (reference) 


20/37 


1.10 (0.46 - 2.61) 




32/41 


1.35 (0.61 - 3.00) 


CRY2 
















rsl 1038689 
















AA 


87/97 


1.0 (reference) 


84/97 


0.99 (0.67 - 1 .47) 




127/127 


1.36 (0.87 - 2.13) 


AG 


52/55 


1.10 (0.66 - 1.81) 


45/57 


0.82 (0.54 - 1 .26) 




68/70 


0.94 (0.65 - 1.37) 


GG 


8/11 


0.86 (0.32 - 2.28) 


13/12 


1.11 (0.49 - 2.51) 

Ptrend' = 0.91 




7/20 


0.42 (0.17 -1.01) 

P trend f = 0.39 


AA 


87/97 


1.0 (reference) 


84/97 


0.99 (0.67 - 1 .47) 




127/127 


1.36 (0.87 - 2.13) 


AG 


52/55 


1.0 (reference) 


45/57 


0.98 (0.58 - 1 .68) 




68/70 


0.82 (0.46 - 1 .47) 


GG 


8/11 


1.0 (reference) 


13/12 


0.44 (0.09 - 2.06) 




7/20 


1 .44 (0.33 - 6.20) 


rs1401417 
















CC 


78/94 


1.0 (reference) 


75/95 


1.01 (0.65 - 1.58) 




1 1 9/1 20 


1.27 (0.84 - 1.91) 


CG 


46/54 


1.10 (0.65-1.86) 


40/56 


0.76 (0.49 - 1.18) 




72/68 


1.13 (0.77 - 1.64) 


GG 


8/8 


1 .36 (0.47-3.94) 


10/13 


0.86 (0.37 - 2.01) 

Ptrend' = 0.61 




4/17 


0.31 (0.10 - 0.94) 

P trend f = 0.78 


CC 


78/94 


1.0 (reference) 


75/95 


1.01 (0.65 - 1.58) 




1 1 9/1 20 


1.27 (0.84 - 1.91) 


CG 


46/54 


1.0 (reference) 


40/56 


0.76 (0.42 - 1 .37) 




72/68 


1.11 (0.65 - 1.90) 


GG 


8/8 


1.0 (reference) 


10/13 


0.22 (0.18 to 3.74) 




4/17 


0.25 (0.05 to 1.37) 


NPAS2 
















rs2305160 
















GG 


75/72 


1.0 (reference) 


58/71 


0.78 (0.48 - 1 .28) 




1 05/90 


1.16 (0.74 - 1.81) 


GA 


61/74 


0.82 (0.50 - 1.33) 


69/83 


0.87 (0.60 - 1 .25) 




94/1 04 


0.93 (0.67 - 1 .29) 


AA 


25/27 


0.89 (0.46 - 1 .72) 


21/28 


0.73 (0.40 - 1 .34) 
P trend f = 0.26 




27/35 


0.87 (0.51 - 1.50) 

Ptrend' = 0.21 


GG 


75/72 


1.0 (reference) 


58/71 


0.78 (0.48 - 1 .28) 




1 05/90 


1.16 (0.74 - 1.81) 


GA 


61/74 


1.0 (reference) 


69/83 


1 .06 (0.65 - 1 .73) 




94/1 04 


1.12 (0.71 - 1.76) 


AA 


25/27 


1.0 (reference) 


21/28 


0.76 (0.33 - 1.72) 




27/35 


0.86 (0.39 - 1 .92) 


ROR-b 
















rs3903529 
















TT 


76/1 00 


1.0 (reference) 


81/89 


1.29 (0.82 - 2.01) 




1 20/1 1 1 


1.41 (0.93 - 2.14) 


TA 


72/58 


1.67 (1.03 - 2.70) 


51/73 


0.74 (0.49 - 1.10) 




81/109 


0.81 (0.58 - 1.13) 
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Table 4 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each polymorph- 
ism among women with less than two consecutive night shifts compared to women with three consecutive night 
shifts or women with four or more consecutive night shifts. (Continued) 



AA 


10/16 


TT 


76/100 


TA 


72/58 


AA 


10/16 


rs7022435 




GG 


89/1 00 


GA 


65/62 


AA 


8/9 


rs7022435 




GG 


89/100 


GA 


65/62 


AA 


8/9 


rs3750420 




cc 


68/87 


CT 


76/72 


TT 


13/15 


cc 


68/87 


CT 


76/72 


TT 


13/15 


AANAT 




rs37601 38 




TT 


38/52 


TG 


86/85 


GG 


40/42 


rr 


38/52 


TG 


86/85 


GG 


40/42 


rs4238989 




GG 


43/60 


GC 


77/82 


CC 


40/35 


GG 


43/60 


GC 


77/82 


cc 


40/35 


AATMD1 A 




rsl 31 13549 




TT 


51/60 


TC 


84/88 


CC 


26/29 


TT 


51/60 


TC 


84/88 


CC 


26/29 


MTNR1B 




rsl 0830963 





0.89 (0.37 - 2.15) 12/19 

0 (reference) 81/89 

.0 (reference) 51/73 

.0 (reference) 1 2/1 9 



.0 (reference) 92/101 
21 (0.75 -1.93) 53/63 
.12 (0.41 - 3.06) 3/18 



.0 (reference) 
.0 (reference) 
.0 (reference) 



.0 (reference) 
.0 (reference) 
.0 (reference) 

.0 (reference) 
.20 (0.72 - 2.0; 
.47 (0.79 - 2.7^ 

.0 (reference) 
.0 (reference) 
.0 (reference) 



.0 (reference) 
19 (0.72-1.96) 
.12 (0.57-2.20) 

.0 (reference) 
.0 (reference) 
0 (reference) 



92/101 

53/63 

3/18 



.0 (reference) 
.0 (reference) 
.0 (reference) 



.0 (reference) 
35(0.84 - 2.18) 68/72 
25 (0.54 - 2.87) 9/25 



68/72 
9/25 



.0 (reference) 42/57 
17 (0.68 - 2.00) 75/95 
.18 (0.63 - 2.21) 



33/32 

42/57 
75/95 
33/32 



0.61 (0.29 - 1.31) 

P ae J = 0.60 

1.29 (0.82 - 2.01) 

0.57 (0.34 - 0.95) 0.03 

0.91 (0.30 - 2.80) 



0.12 0.18 



1 .06 (0.69 - 1 .62) 
0.93 (0.61 - 1.40) 

0.16 (0.05 - 0.55)0.004 0.24 0.23 



Pu 



0.21 



1 .06 (0.69 - 1 .62) 
0.81 (0.48 - 1.37) 

0.17 (0.03 - 0.84)0.03 049 0 28 



1.12 (0.70 - 1.81) 

1.11 (0.76 - 1.62) 

0.40 (0.18 - 0.88) 0.02 

P ae J = 0.58 

1.12 (0.70 - 1.81) 
0.91 (0.56 - 1.48) 
0.34 (0.11 - 1.05) 



0.19 0.22 



0.93 (0.51 - 1.70) 
0.87 (0.61 - 1.24) 
1 .23 (0.72 - 2.08) 
P trend f = 0.03 
0.93 (0.51 - 1.70) 
0.80 (0.51 - 1.25) 
1.28 (0.64 - 2.57) 



50/59 


1.13 (0.64 - 


2.01) 


66/86 


0.87 (0.60 - 


1.27) 


37/39 


1.15 (0.70 - 


1.88) 




P tre J = 0.49 


50/59 


1.13 (0.64 - 


2.01) 


66/86 


0.85 (0.53 - 


1.37) 


37/39 


1 .02 (0.52 - 


2.01) 


64/58 


1 .36 (0.79 - 


2.34) 


67/80 


0.89 (0.61 - 


1.30) 


20/42 


0.50 (0.28 


- 0.9C 




Ptrenc/ - 0.1 1 


64/58 


1 .36 (0.79 - 


2.34) 


67/80 


0.88 (0.55 - 


1.40) 


20/42 


0.55 (0.25 - 


1.23) 



0.11 0.18 



27/15 1.87 (1.00 - 3.61) 0.06 0.24 0.24 

P trend f = 0.28 
120/111 1.41 (0.93 - 2.14) 
81/109 0.66 (0.41 - 1.05) 
27/1 5 2.47 (0.83 - 7.34) 

Ptrend = 0.03 



124/134 1.06(0.71 
82/84 1.05 (0.75 
18/12 



1.57) 
1.49) 

1.55 (0.73 - 3.30) 



0.31 



124/134 1.06 (0.71 - 1.57) 
82/84 
18/12 



93/1 1 2 
1 00/92 
29/21 

93/1 1 2 
1 00/92 
29/21 



69/72 



0.93 (0.58 - 1 .49) 
1 .68 (0.44 - 6.33) 

1.03 (0.66 - 1.61) 
1 .32 (0.95 - 1 .83) 
1.57 (0.87 - 2.83) 

P„ e J = 0.04 

1-03 (0.66 - 1.61) 
1.13 (0.72 - 1.76) 
1.35 (0.51 - 3.58) 



1 9 (0.68 - 2.08) 



108/135 0.91 (0.67 
55/32 



1.24) 



1.82(1.13 - 2.93) 0.01 0.06 0.14 
P trend f = 0.04 
69/72 1.19 (0.68 - 2.08) 
1 08/1 35 0.82 (0.54 - 1 .23) 
55/32 1 .84 (0.94 - 3.62) 

68/84 1 .07 (0.62 - 1 .83) 

120/120 1.18 (0.87 - 1.61) 

43/30 1.62 (1.00 - 2.68) 0.06 0.19 0.22 

P trend f = 0.02 

68/84 1 .07 (0.62 - 1 .83) 

120/120 1.11 (0.73 - 1.68) 

43/30 1 .40 (0.69 - 2.83) 



81/78 1.37 (0.82 - 2.27) 

1 02/96 1 .08 (0.78 - 1 .50) 

47/56 0.82 (0.54 - 1 .26) 

P trend f = 0.95 

81/78 1.37 (0.82 - 2.27) 

102/96 1.11 (0.73 - 1.71) 

47/56 0.79 (0.40 - 1 .59) 
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Table 4 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each polymorph- 
ism among women with less than two consecutive night shifts compared to women with three consecutive night 
shifts or women with four or more consecutive night shifts. (Continued) 



cc 


95/97 


1.0 (reference) 


76/103 


0.74 (0.48 - 


1.13) 


1 1 9/1 29 


0.90 (0.61 - 1 .34) 


CG 


56/69 


0.74 (0.46 - 1.19) 


55/60 


1 .04 (0.69 - 


1.56) 


86/90 


1 .06 (0.76 - 1 .49) 


GG 


7/8 


0.76 (0.25 - 2.32) 


15/13 


1.29 (0.59 - 


2.81) 


19/8 


2.48 (1 .06 - 5.80) 0.04 0 26 0 24 










^trend ~ 0.' 


39 




Ptrenc/ = 0.35 


CC 


95/97 


1.0 (reference) 


76/103 


0.74 (0.48 - 


1.13) 


1 1 9/1 29 


0.90 (0.61 - 1.34) 


CG 


56/69 


1.0 (reference) 


55/60 


1.18 (0.70 - 


2.00) 


86/90 


1 .23 (0.76 - 1 .97) 


GG 


7/8 


1.0 (reference) 


15/13 


1.55 (0.38 - 


6.37) 


19/8 


3.51 (0.73 - 16.96) 



a Number of cases and controls that were successfully genotyped. b OR derived from unconditional logistic regression, including an interaction term between the 
night shift variable and the SNP genotypes and adjusted for age at diagnosis (cases) or interview (controls), selection period, alcohol use, hormone replacement 
therapy for the past 2 years, number of children and history of breast cancer in first-degree family members (mother or sister}. C P interaction was obtained from 
the interaction term between the SNP genotypes and the night shift variable using frequency of the common genotype in the 0 to 2 night shift group as 
reference. Only significant or near-significant values are shown. d The false-positive report probability (FPRP) was calculated as described by Wacholder ef al. [42]. 
FPRP < 0.2 was considered noteworthy. *The Bayesian false discovery probability (BFDP) was calculated as described by Wakefield ef al. [41]. BFDP < 0.8 was 
considered noteworthy. For more details on FPRP and BFDP, see the Methods section. f P tren d values for Cochran-Armitage trend test were calculated from 2xC 
tables using frequency of genotypes in cases and controls from 0 to 2 night shift group as reference. 



carrying the same SNPs and genotypes had a reduced risk 
of breast cancer, indicating a protective role for these 
genetic variants. In addition, women in this exposure 
group carrying at least one variant allele of several SNPs in 
the other core circadian genes such as BMAL2, CSNK1E, 
NPAS2 and PER3 had a noteworthy reduced risk of breast 
cancer. For comparison, only two SNPs in the CRY2 gene 
(rsll038689 and rsl401417) had protective effects in 
women with four or more night shifts, though this effect 
was not noteworthy. 

Many of the associations, such as PER3 (rsl012477), 
CSNK1E (rsl534891) and 0,COr(rslll33373, rs3749474), 
are novel in relation to breast risk but may not be specific 
to breast cancer, because they have also been found to be 
associated with prostate and colorectal cancer risk [47,51]. 
It might be that these SNPs affect expression of the respec- 
tive genes. Current results from the available literature 
on the functionality of other SNPs in these genes are not 
conclusive. For example, it has previously been seen that, in 
healthy nurses and midwives who worked either day or 
rotating night shifts, there were no differences between 
BMAL1, CLOCK, CRY1, CRY2, PERI, PER2 and PER3 gene 
expression levels in peripheral blood leukocytes [52]. In the 
same subjects, however, polymorphisms of the circadian 
genes BMAL1 (rs2279287), CLOCK (rsl801260), PERI 
(rs2735611), PER2 (rs2304672), PER3 (rsl0462020), CRY1 
(rs8192440) and CRY2 (rsl0838527, rsl0838527) were not 
found to be different between shift workers and day-only 
workers [53]. 

In a recent study by Monsees et al. [54], the SNP 
rs2305160 (Ala394Thr) in the NPAS2 gene was found to 
be associated with reduced breast cancer risk in women 
with less than 2 yr of shift work. We were not able to con- 
firm such an association in our study; however, for 
another SNP, rsl7024926 in the NPAS2 gene, we found a 
reduced risk of breast cancer in shift workers who had 



worked three night shifts for at least 5 yr (Table 5). The 
reduced risk was associated with the variant homozygote 
genotype where carriers of the CC variant genotype had 
reduced risk of breast cancer compared to subjects who 
had worked no or few night shifts carrying either the var- 
iant CC genotypes or the common TT genotype (Table 5). 
The NPAS2 rs 17024926 has been found to be associated 
with risk of prostate cancer [47]. This SNP was not geno- 
typed in the study by Monsees et al. [54]; however, the 
reasons why the results of our study and those of Monsees 
et al. differ for Ala394Thr variant are not known. This dif- 
ference could not be due to differences in MAF in the two 
populations, first because the MAF 0.37 for this SNP in 
our population (Additional File 2, Table S2) is similar to 
the MAF 0.34 in the study by Monsees et al., and second 
because we regenotyped this SNP in our subjects and 
again found no associations. There could be several other 
possibilities, such as differences in study design, popula- 
tion differences in incidence of breast cancer, exposure 
metrics and the fact that breast cancer is known to be a 
genetically heterogeneous disease. Nevertheless, both stu- 
dies support a role of circadian genes in breast cancer in 
shift workers, which warrants further confirmation and 
validation in future studies. 

Increased risk of breast cancer in women with four or 
more night shifts was associated with ROR-b SNPs 
rs3903529 and rs3750420 (heterozygotes), whereas women 
with three night shifts and carrying a homozygote variant 
of rs7022435 and rs3750420 had reduced risk. The three 
SNPs are located in intron 1 downstream of the exon 1 in 
the ROR-b transcripts. These SNPs are part of a haplotype 
block that spans the region containing the promoter 
region and the exon 1 (55 kb upstream to 6 kb down- 
stream) of the gene. This haplotype has also been asso- 
ciated with an increased risk of human bipolar syndrome 
[55], but we were not able to find such an association with 
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Table 5 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each 
polymorphism among women with less than two consecutive night shifts for 5 or more years compared to women 
with three consecutive night shifts or women with four or more consecutive night shifts for 5 or more years 



Genes and 
SNPs 



Number of consecutive night shifts nurse has ever worked for 5 or more years 



0 to 2 consecutive 
nights 

Ca a / OR (95% Cl) b Ca a / 
Co a Co a 



3 consecutive nights 



4 or more consecutive nights 



OR (95% Cl) b P in 



FPRP" 
BFDP e 



Ca a / OR (95% Cl) b P inter c 
Co a 



FPRP" 
BFDP e 



ARNTL 1/ 
BMAL1 

rs2290035 

TT 

TA 

AA 



TT 
TA 
AA 

rs2278749 

CC 

CT 



TT 



CC 
CT 
TT 

rs969485 

AA 

AG 

GG 



AA 
AG 
GG 

ARNTL2/ 
BMAL2 

rs2306074 

TT 

TC 

CC 



TT 
TC 
CC 

CLOCK 

rs3749474 
CC 



106/121 
149/161 



50/61 



106/121 
149/161 
50/61 

212/231 
86/87 

7/19 



212/231 

86/87 

7/19 



1.0 (reference) 

1 .03 (0.72 - 
1.47) 

0.95 (0.59 - 
1.52) 

1.0 (reference) 
1.0 (reference) 
1.0 (reference) 

1.0 (reference) 

1.19 (0.83 - 
1.73) 

0.40 (0.16 - 
0.98) 

1.0 (reference) 
1.0 (reference) 
1.0 (reference) 



161/189 1 .0 (reference) 
1 05/1 17 1 .05 (0.74 - 
1.49) 

21/20 1.07 (0.54 - 
2.09) 

161/189 1 .0 (reference) 
1 05/1 17 1 .0 (reference) 
21/20 1 .0 (reference) 



143/144 1 .0 (reference) 

1 34/1 63 0.84 (0.60 - 
1.17) 



37/39 



143/144 
1 34/1 63 
37/39 



1.04 (0.61 - 
1.76) 

1 .0 (reference) 
1.0 (reference) 
1.0 (reference) 



40/46 
66/71 



0.98 (0.58 - 1 .65) 
1.04 (0.71 - 1.51) 



1 6/27 0.69 (0.36 - 1 .34) 

P trend f = 0.60 
40/46 0.98 (0.58 - 1 .65) 
66/71 0.99 (0.65 - 1.50) 
1 6/27 0.75 (0.35 - 1 .60) 

97/92 1.19 (0.84 - 1.70) 

22/45 0.51 (0.29 - 0.87) 0.01 0.07 0.15 

1/5 0.26 (0.03 - 2.38) 

P trend f = 0.07 
97/92 1.19 (0.84 - 1.70) 
22/45 0.47 (0.25 - 0.88) 0.02 0.12 0.18 

1/5 1.33 (0.10 - 17.78) 

51/74 0.83 (0.54 - 1.28) 
58/55 1 .23 (0.82 - 1 .85) 

8/8 1 .40 (0.50 - 3.94) 

P ae „J = 0.49 

51/74 0.83 (0.54 - 1.28) 
58/55 1 .50 (0.72 - 1 .85) 
8/8 1 .35 (0. 37 - 4.88) 



58/63 0.91 (0.59 - 1.43) 

60/68 0.97 (0.66 - 1 .42) 

6/21 0.30 (0.17 - 0.75) 0.01 0.18 0.21 

P trend f = 0.06 

58/63 0.91 (0.59 - 1.43) 

60/68 1 .06 (0.68 - 1 .64) 

6/21 0.24 (0.08 - 0.77) 0.02 0.31 0 25 



26/26 
46/50 

36/21 



26/26 
46/50 
36/21 

66/58 
33/32 

5/6 



66/58 
33/32 
5/6 



0.99 (0.53 - 1 .88) 
1.02 (0.66 - 1.57) 

1.91 (1.08 - 3.37) 0.03 0.12 0.19 

P trend f = 0.08 
0.99 (0.53 - 1.88) 
1.10 (0.64 - 1.63) 

1.95(0.97 - 3.90) 0.06 0.25 0.24 



1 .26 (0.83 ■ 
1.10 (0.64 ■ 



.91) 
.84) 



0.84 (0.25 - 2.82) 

Pvend ~ 0-52 

1.26 (0.83 - 1.91) 
0.99 (0.55 - 1.79) 
2.68 (0.48 - 15.01) 



41/50 0.92(0.57 - 1.49) 

50/35 1.64(1.03 - 2.61) 0.04 0.12 0.18 

9/7 1 .50 (0.53 - 4.22) 

P trend f = 0.05 

41/50 0.92 (0.57 to 1.49) 

50/35 1.63(0.96 - 2.76) 0.07 0.21 0.23 

9/7 1 .60 (0.38 - 6.74) 



51/53 0.97 (0.61 - 1.53) 

50/40 1.36(0.86 - 2.13) 

9/9 1 .03 (0.39 - 2.72) 

P vend = 0-5 1 

51/53 0.97 (0.61 - 1.53) 

50/40 1.37(0.83 - 2.26) 

9/9 0.99(0.32 - 3.11) 



123/131 1 .0 (reference) 52/55 0.97 (Oi 



1.37) 



48/42 1 .08 (0.73 - 1 .59) 
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Table 5 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each polymorph- 
ism among women with less than two consecutive night shifts for 5 or more years compared to women with three 
consecutive night shifts or women with four or more consecutive night shifts for 5 or more years (Continued) 



CT 


148/152 


1 .04 (0.74 - 
1.47) 


56/63 


0.87 (0.46 - 1 .63) 




47/45 


1 .84 (0.92 - 3.69) 


TT 


39/66 


0.60 (0.37 - 
0.97) 


16/29 


0.24 (0.05 - 1.18) 0.07 

P trend f = 0.24 


0.56 0.29 


13/13 


3.15 (0.62 - 16.14) 
P trend f = 0.63 


CC 


123/131 


1.0 (reference) 


52/55 


0.97 (0.68 - 1 .37) 




48/42 


1.08 (0.73 - 1.59) 


CT 


148/152 


1.0 (reference) 


56/63 


0.64 (0.30 - 1 .36) 




47/45 


1.44 (0.63 - 3.31) 


TT 


39/66 


1.0 (reference) 


16/29 


0.30 (0.05 - 1 .77) 




13/13 


2.12 (0.30 - 15.15) 


CSNK1E 
















rs1 534891 
















GG 


241/258 


1.0 (reference) 


98/1 1 2 


0.93 (0.66 - 1 .309) 




82/76 


1.12 80.77 - 1.62) 


GA 


62/77 


0.83 (0.57 - 
1.23) 


17/35 


0.55 (0.30 - 1.01) 0.05 


0.21 0.23 


25/25 


1.09 (0.61 - 1.95) 


A A 




0.51 (0.16 - 
1.60) 


7/1 
Z/ 1 


i nn in 1 7 93 en 1 ! 

Ptrenc/ =0.13 




i n 
i /z 


ftrend' = 0.66 


GG 


241/258 


1.0 (reference) 


98/1 1 2 


0.93 (0.66 - 1 .309) 




82/76 


1.12 80.77 - 1.62) 


GA 


62/77 


1.0 (reference) 


17/35 


0.63 (0.31 - 1.28) 




25/25 


1.11 (0.56 - 2.21) 


AA 


5/9 


1.0 (reference) 


2/1 


2.58 (0.14 - 91.16) 




1/2 


2.84 (0.08 - 96.47) 


rs5757037 
















CC 


126/154 


1.0 (reference) 


49/57 


1.10 (0.68 - 1.76) 




41/38 


1.28 (0.75 - 2.70) 


CT 


137/140 


1.11 (0.79 - 
1.58) 


51/77 


0.67 (0.45 - 0.99) 0.05 


0.13 0.18 


49/45 


1.14 (0.74 - 1.77) 


TT 


43/43 


1 .04 (0.63 - 
1.71) 


14/16 


0.90 (0.43 - 1 .90) 




16/16 


1.00 (0.49 - 2.05) 








Ptrend' = 0.57 






Ptrend' = 0.22 


CC 


126/154 


1.0 (reference) 


49/57 


1.10 (0.68 - 1.76) 




41/38 


1.28 (0.75 - 2.70) 


CT 


137/140 


1.0 (reference) 


51/77 


0.65 (0.42 - 1 .02) 




49/45 


1.08 (0.66 - 1.75) 


TT 


43/43 


1.04 (0.63- 1.71) 


14/16 


0.89 (0.38 - 2.10) 




16/16 


0.97 (0.42 - 2.23) 


NPAS2 
















rs2305160 
















GG 


1 36/1 34 


1.0 (reference) 


52/57 


0.86 (0.54 - 1 .36) 




50/42 


1.15 (0.70 - 1.88) 


GA 


130/151 


0.84 (0.59 - 
1.18) 


51/69 


0.84 (0.56 - 1 .25) 




43/41 


1.10 (0.70 - 1.74) 


AA 


41/56 


0.72 (0.45 - 
1.18) 


18/21 


0.87 (0.45 - 1 .69) 

P trend = 0-2 1 




14/13 


1.22 (0.56 - 2.64) 
P trend f = 0.79 


rs2305160 
















GG 


1 36/1 34 


1.0 (reference) 


52/57 


0.86 (0.54 - 1 .36) 




50/42 


1.15 (0.70 - 1.88) 


GA 


130/151 


1.0 (reference) 


51/69 


0.95 (0.61 - 1.50) 




43/41 


1.19 (0.72 - 1.96) 


AA 


41/56 


1.0 (reference) 


18/21 


1.10 (0.50 - 2.40) 




14/13 


1.43 (0.58 - 3.52) 


rs 17024926 
















TT 


1 1 0/1 33 


1.0 (reference) 


60/66 


1 .22 (0.78 - 1 .92) 




36/45 


1.03 (0.61 - 1.75) 


TC 


148/178 


1 .08 (0.76 - 
1.52) 


55/64 


0.98 (0.66 - 1 .47) 




54/46 


1 .23 (0.80 - 1 .89) 


CC 


47/38 


1.51 (0.90 - 
2.52) 


6/21 


0.33 (0.13 - 0.84) 
0.020 

P tren d f = 0.28 


0.24 0.23 


17/10 


1.89 (0.85 - 4.23) 

P tre nd f = 0.05 


TT 


1 1 0/1 33 


1.0 (reference) 


60/66 


1 .22 (0.78 - 1 .92) 




36/45 


1.03 (0.61 - 1.75) 


TC 


148/178 


1.0 (reference) 


55/64 


1.01 (0.65 - 1.57) 




54/46 


1.28 (0.81 - 2.04) 


CC 


47/38 


1.0 (reference) 


6/21 


0.24 (0.09 - 0.67) 


0.19 0.22 


17/10 


1.49 (0.56 - 3.95) 



0.006 



PER2 
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Table 5 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each polymorph- 
ism among women with less than two consecutive night shifts for 5 or more years compared to women with three 
consecutive night shifts or women with four or more consecutive night shifts for 5 or more years (Continued) 



rs 11695472 

AA 

AC 

CC 



AA 
AC 
CC 
PER3 

rsl 01 2477 

GG 

GC 



CC 



GG 
GC 

CC 

ROR-b 

rs3903529 

TT 

TA 

AA 



TT 
TA 
AA 

rs7022435 

GG 

GA 

AA 



GG 
GA 
AA 

rs3750420 

CC 

CT 

TT 



CC 

CT 

rr 



1 73/1 87 1 .0 (reference) 
117/143 0.90 (0.65 - 
1.25) 



22/21 



1 73/1 87 
117/143 
22/21 



203/237 
95/1 00 

9/7 



203/237 
95/100 



9/7 



1 62/1 83 
123/131 



22/27 



1 62/1 83 
123/131 
22/27 

1 78/1 94 
115/128 

15/19 



1 78/1 94 
115/128 
15/19 

136/166 
136/143 

28/31 



1 36/1 66 
136/143 
28/31 



.24 (0.64 ■ 
.39) 



.0 (reference) 
.0 (reference) 
.0 (reference) 



.0 (reference) 

.10 (0.78 - 
■56) 

.69 (0.58 - 
.93) 

.0 (reference) 
.0 (reference) 



.0 (reference) 

.08 (0.77 - 
51) 

.98 (0.53 - 
.82) 

.0 (reference) 
.0 (reference) 
.0 (reference) 

.0 (reference) 

.00 (0.72 - 
.40) 

.93 (0.45 - 
.90) 

.0 (reference) 
.0 (reference) 
.0 (reference) 

.0 (reference) 

.20 (0.86 - 
.68) 

.18 (0.67 - 



.0 (reference) 
.0 (reference) 
.0 (reference) 



68/87 
38/54 



16/7 



68/87 
38/54 
16/7 



88/97 
22/51 



3/1 



88/97 
22/51 



.0 (reference) 3/1 



64/70 
41/62 



0.87 (0.59 - 1 .29) 
0.74 (0.47 - 1.16) 

2.69 (1.08 - 6.73) 0.03 0 28 0 25 

P ae J = 0.93 
0.87 (0.59 - 1 .29) 
0.83 (0.50 - 1 .39) 
1.91 (0.59 - 6.22) 



1.12 (0.77 - 1.58) 

0.42 (0.24 - 0.72) 
0.002 

3.42 (0.35 - 33.40) 



Ptrend = 0.15 

1.12 (0.77 - 1.58) 

0.38 (0.18 - 0.66) 
0.001 

3.00 (0.12 - 77.00) 



1 .09 (0.71 - 1 .66) 
0.75 (0.49 - 1.15) 



0.02 0.08 



0.01 0.06 



1 3/1 8 0.76 (0.36 - 1 .59) 

P trend f = 0.26 
64/70 1.09 (0.71 - 1.66) 
41/62 0.72(0.44 -1.16) 
13/18 0.78 (0.30 - 2.05) 

72/85 0.97 (0.65 - 1 .43) 
44/49 1 .02 (0.65 - 1 .60) 

3/12 0.25 (0.07 - 0.89) 0.03 0 410 27 

Ptrend' = 0.25 

72/85 0.97 (0.65 - 1 .43) 
44/49 1 .03 (0.62 - 1 .70) 
3/1 2 0.29 (0.06 - 1 .33) 

52/66 0.99 (0.63 - 1 .55) 
60/58 1 .24 (0.83 - 1 .86) 

6/18 0.37 (0.14- 0.96) 0.04 0.310 25 

P, rend f = 0.82 
52/66 0.99 (0.63 - 1 .55) 
60/58 1.13 (0.72 - 1.77) 
6/18 0.33 (0.11 - 1.00) 0.05 0.39 0 27 



63/48 
36/44 

9/8 



63/48 
36/44 
9/8 



3/2 



35/24 
3/2 



51/47 
40/47 

14/5 



51/47 
40/47 
14/5 

55/56 
41/32 

11/8 



55/56 
41/32 
11/8 



1.37 (0.88 - 2.14) 
0.87 (0.54 - 1 .39) 

1.32 (0.48 - 3.65) 

Ptrend' = 0.79 

1.37 (0.88 - 2.14) 
0.93 (0.55 - 1.55) 
0.83 (0.22 - 3.14) 



66/73 1.02(0.69 - 1.51) 
35/24 1.56 (0.91 - 2.68) 



1.90 (0.31 



.58) 



P,rend = 0.08 

66/73 1.02(0.69 -1.51) 
1.46 (0.79 - 2.72) 



1.30 (0.06 - 28.93) 



1 .20 (0.75 ■ 
0.92 (0.59 ■ 



1 .93) 
1 .46) 



2.77 (0.98 - 7.86) 
0.056 

Ptrend' = 0.21 

1.20 (0.75 - 1.93) 
0.94 (0.57 - 1.55) 
3.07 (0.84 - 11.18) 

1 .05 (0.67 - 1 .64) 
1.38 (0.84 - 2.26) 

1.34 (0.52 - 3.41) 

Ptrend' =0.15 

1 .05 (0.67 - 1 .64) 
1.34 (0.78 - 2.92) 
2.18 (0.55 - 8.61) 



0.38 0.27 



41/47 1.02 (0.61 - 1.68) 

48/35 1.60(1.00 - 2.56) 0.05 0.15 0.20 

17/12 1.52 (0.71 - 3.28) 



41/47 
48/35 
17/12 



P, rend ' = 0.03 
1.02 (0.61 - 1.68) 
1.47 (0.88 - 2.45) 
1.39 (0.54 - 3.60) 
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Table 5 Odds ratios and 95% confidence intervals comparing odds of breast cancer associated with each polymorph- 
ism among women with less than two consecutive night shifts for 5 or more years compared to women with three 
consecutive night shifts or women with four or more consecutive night shifts for 5 or more years (Continued) 



MTNR1A 


















rs1 31 13549 


















TT 


1 1 3/1 09 


1.0 (reference) 


48/54 


0.94 (0.57 - 


1.53) 


35/33 


1 .09 (0.62 - 1 .93) 




TC 


142/155 


0.94 (0.66 - 


59/60 


1 .08 (0.72 - 


1.61) 


52/49 


1.08 (0.71 - 1.65) 








1.35) 














CC 


55/75 


0.75 (0.48 - 
1.18) 


16/34 


0.47 (0.25 


- 0.91) 0.02 0.15 0.20 


22/18 


1.24 (0.65 - 2.39) 












Ptrenc/ = 0.C 


8 




Ptrenc/ = 0.71 




TT 


1 1 3/1 09 


1 .0 (reference) 


48/54 


0.94 (0.57 - 


1.53) 


35/33 


1.09 (0.62 - 1.93) 




TC 


142/155 


1.0 (reference) 


59/60 


1 .03 (0.66 - 


1.61) 


52/49 


1.09 (0.68 - 1.75) 




CC 


55/75 


1.0 (reference) 


16/34 


0.57 (0.27 - 


0.19) 


22/18 


1.28 (0.61 - 2.68) 




MTNR1B 


















rs1 0830963 


















CC 


1 76/1 95 


1.0 (reference) 


63/84 


0.83 (0.58 - 


1.23) 


51/50 


1.01 (0.64 - 1.61) 




CG 


105/123 


0.86 (0.61 - 
1.21) 


46/51 


1 .06 (0.69 - 


1.64) 


46/45 


1.13 (0.73 - 1.75) 




GG 


21/21 


1.05 (0.55 - 


11/7 


1.59 (0.59 - 


4.25) 


9/1 


9.12 (1.14 - 73.09) 


0.59 0.29 






2.02) 










0.04 












Ptrenc/ = 0.76 




Ptrenc/ = 0.10 




CC 


176/195 


1.0 (reference) 


63/84 


0.83 (0.58 - 


1.23) 


51/50 


1.01 (0.64 - 1.61) 




CG 


105/123 


1.0 (reference) 


46/51 


1.17 (0.71 - 


1.93) 


46/45 


1.26 (0.76 - 2.08) 




GG 


21/21 


1.0 (reference) 


11/7 


1 .44 (0.43 - 


4.86) 


9/1 


9.45 (0.97 - 91.82) 





a Number of cases and controls that were successfully genotyped. b OR derived from unconditional logistic regression, including an interaction term between the 
night shift variable and the SNP genotypes, and adjusted for age at diagnosis (cases) or interview (controls), selection period, alcohol use, hormone replacement 
therapy for the past 2 years, number of children and history of breast cancer in first-degree family members (mother or sister). C P interaction was obtained from 
the interaction term between the SNP genotypes and the night shift variable using frequency of the common genotype in the 0 to 2 night shift group as 
reference. Only significant or near-significant values are shown. d The false-positive report probability (FPRP) was calculated as described by Wacholder et al. [42]. 
FPRP < 0.2 was considered noteworthy. *The Bayesian false discovery probability (BFDP) was calculated as described by Wakefield ef al. [41). BFDP < 0.8 was 
considered noteworthy. For more details on FPRP and BFDP, see the Methods section. f P tr end values for Cochran-Armitage trend test were calculated from 2xC 
tables using frequency of genotypes in cases and controls from 0 to 2 night shift group as reference. 



breast cancer risk (data not shown). If these SNPs, or 
others in linkage disequilibrium with them, may represent 
susceptibility variants for breast cancer, it is possible they 
do so by altering the gene expression of the ROR-b gene. 
The ROR proteins are involved in regulation of circadian 
rhythms by activating transcription of the core circadian 
BMAL1 gene [56]. 

Two polymorphisms in the CRY2 gene (rsl 1038689 and 
rsl401417) had protective effects in women with four or 
more night shifts. The biological functions of these SNPs 
are not known. These polymorphisms have been found to 
be associated with postmenopausal breast cancer risk, 
with significant effect modification by menopausal status 
where this association was evident only in women with 
ER- and/or PR-negative breast tumors but not ER- and/or 
PR-positive tumors [57]. This may indicate that these 
polymorphisms may influence breast cancer risk through 
estrogen and progesterone signaling pathways. ERs and 
PRs are important for reproductive performance in 
women. Kovanen et al. [46] showed that women with 
BMAL1 rs2278749 TT genotypes had more pregnancies 
and miscarriages, whereas NPAS2 rsl 1673746 T carriers 



had fewer miscarriages, indicating a possible role of these 
polymorphisms affecting age at first birth and also number 
of children, both of which are established risk factors for 
breast cancer [58]. 

In the genes controlling melatonin biosynthesis (AANA T) 
or melatonin receptors (MTNR1A and MTNR1B), several 
SNPs were associated with increased risk of breast cancer. 
In the main effects analysis, the variant genotypes of both 
of the AANAT polymorphisms were associated with 
increased breast cancer risk. However, the risk associated 
with the rs3760138 GG genotype was 40% higher in the 
most exposed group with four or more night shifts (1.82 vs 
1.42). These two SNPs are located in the regulatory regions 
of the AANA T gene which could potentially affect the 
expression of the AANAT gene. The AANAT gene product 
is an acetyltransferase that regulates melatonin biosynthesis 
during night and day periods by pineal gland. Polymorph- 
isms in this gene have also been associated with sleep disor- 
ders, but there are no data on the association between these 
polymorphisms and breast cancer risk. The levels of mela- 
tonin may also be affected by variants in the genes that 
encode melatonin receptors 1A and IB. The MTNR1B 
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rsl0830963 was associated with increased risk in the high- 
level exposed group of women carrying two variant alleles, 
without being noteworthy by FPRP and BFDP. The variant 
alleles of MTNR1A rsll3113549 polymorphisms were asso- 
ciated with decreased risk of breast cancer in the women 
with three consecutive night shifts. Night shift workers 
have been shown to have substantially reduced melatonin 
levels during night work and daytime sleep, and levels 
remain low even when a night shift worker sleeps at night 
[22]. 

Strengths and limitations of the study 

Our results indicate that some SNPs of the core circa- 
dian or melatonin signaling pathways, in combination 
with night work, may affect breast cancer risk. These 
findings are consistent with the results of recent studies, 
including our previous study of the same nurses, sug- 
gesting a higher risk with longer duration of intense 
night shift work per week [7,9,10]. The reported associa- 
tions are biologically plausible, as the polymorphisms in 
the genes investigated could influence gene expression, 
protein function and protein-protein interactions, all of 
which may affect several known risk factors for breast 
cancer. However, this is a medium-sized, exploratory 
study highlighting the need for investigation of possible 
combinatory effects of genetic polymorphisms and night 
shift work in breast cancer. The results should be inter- 
preted with caution, and larger studies are needed to 
confirm the associations found here. 

The study may have a number of limitations. A larger 
reference group without any night shift work is needed. 
Our decision to use women with maximum of two conse- 
cutive night shifts as the reference group was made to 
ensure a sufficient number of subjects for each genotype. 
Furthermore, the metrics of exposure to night work used 
in our study are different from those used in other studies; 
therefore, the results may not be comparable. We have 
selected a limited number of polymorphisms, and, to 
cover all SNPs in the core circadian and melatonin signal- 
ing pathways, a much larger study (that is, GWAS) is war- 
ranted. The selection of genes and polymorphisms differs 
from other studies as well and may have affected the out- 
comes. Finally, although we set very low prior probabilities 
for multiple testing, it is still possible that some associa- 
tions may not be true positive or negative findings. 

Conclusion 

We have found a number of polymorphisms in a num- 
ber of core circadian or melatonin biosynthesis and/or 
bioavailability pathway genes that could be associated 
with breast cancer risk in shift workers. The associations 
seem biologically plausible because many of the poly- 
morphisms may alter gene expression, mRNA stability 
or protein function. Future studies investigating the 



genetic and epigenetic functionality of these SNPs, as 
well as epigenetic changes in the shift workers, are war- 
ranted to understand exact biological mechanisms. 
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